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DPMJET is a Monte Carlo model for sampling of hadron-hadron, hadron-nucleus, nucleus- 
nucleus and neutrino-nucleus interactions at accelerator and Cosmic Ray energies according to the 
two-component Dual Parton Model 

Here we describe new features in version DPMJET-II.5: Implementation of new DPM diagrams 
for an improved description of baryon stopping in nuclear collisions and improvements in the cal- 
culation of Glauber cross sections. The new diagrams allow two quite different extrapolations of 
the model to the highest Cosmic Ray energies. The new version of the model is compared to 
■ experimental data on hadron-hadron, hadron-nucleus and nucleus-nucleus collisions. 
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I. INTRODUCTION 



In the present paper we desribe new features in the physics of DPMJET-II.5 and compare important features of 
the model to experimental data. The code manual for DPMJET-II.5 will be given in a separate paper jl]. 
' The DPMJET-II.2 event generator was described in detail before [|[jD. The main topic in switching to DPMJET- 
^ ' II.3/II.4 H was the extension of the model up to energies of = 2000 TeV. The two versions DPMJET-II.3 and 
fSj DPMJET-II.4 did differ mainly by the versions of the Lund code JETSET they use. DPMJET-II.3 did use a JETSET 
I : version transformed into DOUBLE PRECISION. DPMJET-II.4 uses the 1997 version PYTHIA-6.1 now in DOUBLE 
^ . PRECISION in which JETSET is contained. DPMJET-II.5 continues to use PYTHIA-6.1. The extension to higher 
energies was done in DPMJET-II.3/4 by calculating the minijet component of the model using as default the GRV- 
^ ■ L094 parton distributions [||. This is replaced in DPMJET-IL5 by the more recent version GRV-LO-98 0. These 
f| ' new parton distributions describe the data measured in the last years at the HERA Collider. 

Qh! DPMJET-II is used in Collaboration with Battistoni, Carboni and Forti |^ for the simulation of the Cosmic Ray 
• cascade in the HEMAS-DPM code system. A hadron production model to be used for simulations at accelerator and 
\ Cosmic Ray energies should use all possible information from fixed target experiments and collider experiments at 
■ accelerators. There might be large differences on what is considered important: Best studied at accelerators is the 
^ ' central region in hadron-hadron collisions. For studying the Cosmic Ray cascade, the main interest is in the forward 
• i-H , fragmentation region of hadron-nucleus and nucleus-nucleus collisions. DPMJET simulates hadron production in the 
framework of the Dual Parton Model with emphasis as well in the central as in the fragmentation region. 
^ , Soft multiparticle production characterizing hadronic interactions at supercollider or Cosmic Ray energies cannot be 
■ - - ' understood purely within theoretical approaches provided by perturbative QCD. The nonperturbative soft component 
of hadron production, which is responsible for all of hadron production at low energies is still acting at higher energies 
together with the hard component described by perturbative QCD. 

Using basic ideas of the dual topological unitarization scheme j|,|l^ the Dual Parton Model (DPM) ||ll| (a published 
review stressing mainly non-Monte Carlo applications is given in Ref. [Q, a nice review, also with emphasis to Monte 
Carlo implementations of the DPM-model and to photon induced reactions can be found in the Ph.D Thesis of Ralph 
Engel ) has been very successfully describing soft hadronic processes together with the gradual transition to hard 
collisions with rising energy. 

Observations like rapidity plateaus and average transverse momenta rising with energy, KNO scaling violation, 
transverse momentum- multiplicity correlations and minijets pointed out, that soft and hard processes are closely 
related. These properties were understood within the two-component Dual Parton Model p^-pO|. 

The hard component is introduced applying lowest order of perturbative hard constituent scattering [|l| . Single 
diffraction dissociation is represented by a triple-pomeron exchange (high mass single diffraction) and a low mass 
component. 

The Dual Parton model provides a framework not only for the study of hadron-hadron interactions, but also for 
the description of particle production in hadron-nucleus and nucleus-nucleus collisions at high energies. Within 
this model the high energy projectile undergoes a multiple scattering as formulated in Glaubers approach; particle 
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production is realized by the fragmentation of colorless parton-parton chains constructed from the quark content of 
the interacting hadrons and nuclei. 

The first successfull applications of the Monte Carlo version of the dual parton model (DPM) to hadron-nucleus 
p2| , ^ and nucleus-nucleus |2^-|2^] collisions also demonstrated that the cascade of created secondaries in the target 
(and projectile) nuclei contributes significantly to particle production in the target (and projectile) fragmentation re- 
gions. DPM JET uses the concept of a formation zone |2^j2^ suppressing in a natural way the cascading of high-energy 
secondaries. The Monte Carlo model includes intranuclear cascade processes of the created secondaries combined with 
the formation time concept since the first version of DPMJET-II. The FZIC (formation zone intranuclear cascade) was 
also introduced for the projectile nucleus and in addition the nuclear evaporation and fragmentation of the residual 
nucleus was introduced into DPMJET-II [|9|j3^. 

In the following we briefly sketch new features of the model and mention the most important ingredients; for a more 
detailed description of the model as applied in the code we refer to Refs. [|J|j2|,||,|l|-|3|j2|@|| The application 
of DPMJET to neutrino-nucleus interactions is desribed in Q . In Section II we will describe new features of the Dual 
Parton Model as used in DPMJET-II. 5. The most important of these new features are new diagrams contributing to 
baryon stopping and a better calculation of Glauber cross sections. In Section III we compare the model to data and 
in Section IV we discuss the properties of the model at the highest energies. In both these Sections we present only 
material obtained with DPMJET-II. 5. It is not to be expected that in the energy region where experimental data 
are available the model has changed so strongly, that comparisons given in former papers with previous versions of 
the model would have become invalid. 

What is new in DPMJET II.5? 

Here we only summarize the new features of DPMJET-II. 5 against DPMJET-II. 3/4 |^ . which are described in 
more detail in later Sections. 

Modifications of the Glauber model calculation of nuclear cross sections. For some light nuclei the 
Woods-Saxon nuclear densities are replaced by parametrizations in better agreement to data and measured nuclear 
radii are used Q instead of the nuclear radii according to a parametrization used before. The new option XSECNUC 
is added to calculate a table of total, elastic, quasi-elastic and production cross sections using a modified version of 
the routine XSGLAU adopted from DTUNUC-II Q,^. These changes lead to an improved agreement to measured 
nuclear cross sections , especially the p-Air cross sections measured in Cosmic Ray experiments. 

Elimination of some old features of the code. Originally, DPMJET was written in single precision and 
compiler options were used on 32 bit word computers to transform the code to DOUBLE PRECISION. Now DPMJET 
is written explicitely in DOUBLE PRECISION. A number of features of the code was changed to get it running on 
different UNIX systems, especially HP-unix, DEC-alphas, IBM-unix, Sun-workstations and LINUX. 

Old features of the code were removed, in particular: the optional fragmentation of strings according to the 
BAMJET code Only the PYTHIA (JETSET)string fragmentation remains in the present DPMJET, we use 

PYTIIIA-6.1 in the DOUBLE PRECISION version. However, the internal particle codes used are still the ones from 
BAMJET, The output particle codes are as before the ones according to the PDG convention. Some of the older 
parton structure functions are dropped from the code. The chain fusion option, which was no longer up to date, has 
been removed. 

Single diffraction also implemented in nucleus nucleus collisions. The present version of the model 
DPMJET-II.5 allows to include or exclude single diffractive events also in the case of nucleus-nucleus collisions or to 
sample only single diffractive events. The diffractive cross section in hadron-nucleus and nucleus-nucleus collisions is 
calculated by the Monte Carlo method. 

Implementation of new DPM diagrams for an improved description of baryon stopping in hadron- 
nucleus and nucleus— nucleus collisions. A new striking feature of hadron production in nuclear collisions is the 
large stopping of the participating nucleons (compared to hadron-hadron collisions) in hadron-nucleus and nucleus- 
nucleus collisions [^,^. The popcorn mechanism implemented in models with independent string fragmentation 
like the DPM since more than 15 years is not sufficient to explain this enhanced baryon stopping. New DPM- 
diagrams were proposed by Capella and Kopcliovich and investigated in detail by Capella and Collaborators 
and by Casado DPMJET-II.5 implements these diagrams and obtains an improved agreement to net- 
baryon distributions in nuclear collisions. 

Implementation of the new diagrams allows for two different extrapolations of the model to Cosmic 
Ray energies beyond the accelerator energy range. The diagram investigated by Capella and Collaborators 
pl| , ^ uses sea quarks, which in nuclear collisions according to the Glauber model are needed to implement the 
multiple collisions in the nucleus. In the Dual Parton Model in the collider energy range there are further multiple 
collisions, again implemented using sea quarks at the ends of multiple chains, which occur due to the unitarization 
procedure. 

We get at high energy enhanced baryon stopping even in proton-proton collisions and a strong modification of the 
model if also these sea quarks are used in the new diagrams The two models obtained with or without this effect 
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result in quite different model predictions at the highest Cosmic Ray energies. 

Improvements in the sampling and analysis of central collisions. Older versions of DPMJET were also 
applied to central heavy ion collisions at CERN-SPS, RHIC and LHC energies. The previous version of DPMJET 
DPMJET-II.3/4 |], was mainly implemented with the aim to sample minimum bias hadron-nucleus and nucleus- 
nucleus collisions for the simulation of the Cosmic Ray cascade. Heavy targets or projectiles like Pb or Au and central 
collisions are rather unimportant for this purpose. 

Many new, previously unknown features of multiparticle production have been found by the experimental heavy 
ion program at the CERN-SPS. It is certainly useful to use these results as a check even if the aim would be just 
to simulate Cosmic Ray events. Important features of these results concern in particular strangeness production and 
the baryon stopping discussed already above in central nucleus-nucleus collisions. 

DPMJET-II.5 does not desribe all features of the heavy ion experiments, nevertheless, it has been tuned to work 
for central collisions and to describe already quite a number of features of the CERN-SPS heavy ion data. 



II. BASIC PHYSICAL CONCEPTS OF THE DUAL PARTON MODEL 



A. The two— component Dual Parton Model for hadron— hadron collisions 



The two component Dual Parton model for hadron-hadron collisions as applied in DPMJET has been described in 
detail in [|5|,p|,p^ pO| . We will not repeat this here. 

The only topic which has to be discussed is the fit of the model parameters using the new parton structure function 
GRV-LO-98 §. 
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FIG. 1. Total, inelastic and elastic pp and pp cross sections as function of the center of mass energy ^/s. The model results 
obtained using the GRV-L098 parton distributions Q are compared to the Donnachie-Landshoff fit for the total cross section 
Q and to data 
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FIG. 2. Single diffractive -pp and pp cross sections as function of the center of mass ener; 
using the GRV-L098 parton distributions |M] are compared to data for M^/s = ^ < 0.05 
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IS. The model results obtained 
to data. 



We fit the Pomeron parameters with the energy dependent cut-off p(](| for minijets: 
Vy-thr = 2.5 + 0.12[logio(Vi/y^o)]'' [GeV/c], = 50GeV. 



(2.1) 



To describe the high energy particle production we have to determine the free parameters of the model, i.e. the 
proton- Pomeron coupling constant the effective soft Pomeron intercept a(0), the slope of the Pomeron trajectory 
a', the slope parameters h and and the excitation coupling constant A. This has been done by a global fit to all 
available data of total, elastic, inelastic, and single diffractive cross sections in the energy range from ISR to collider 
experiments as well as to the data on the elastic slopes in this energy range. Since there are some differences in the 
hard parton distribution functions at small x values resulting in different hard input cross sections we have to perform 
separate fits for each set of parton distribution functions. 

We get good fits using the GRV-LO-98 PDF. In Table | we give the parameters obtained in the fit. The values 
obtained for a(0) demonstrate again that the fits result in a supercritical Pomerom. In Fig. |l| we plot the fitted 
cross sections together with the data. Furthermore we compare the total cross sections obtained with the popular 
Donnachie-Landshoff fit [Q . In Fig. || we plot the fitted single diffractive cross sections together with data [p4|-p8[ . 

TABLE I. DTUJET99 model parameters obtained with an energy dependent p±thr ■ 
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The single diffractive cross section becomes energy independent at high energies in DTUJET and DPMJET. This 
foUows from the triple Pomeron cross section chosen in |15[| without the energy rise according to the supercritical 
Pomeron. (As explained in detail elsewhere DPMJET follows in the unitarization procedure DTUJET.) Such an 
approximation is not necessary in the Dual Parton Model. It has been shown by Engel within the PHOJET Dual 
Parton Model |^9| that the proper inclusion of the supercritical triple Pomeron cross section as well as higher order 
terms like the double pomeron graph into the unitarization procedure also results in a single diffractive cross section 
becoming rather energy independent at high energy in agreement to the data. Similar results were obtained by 
Gotsman, Levin and Maor |Q . On a phenomenological basis this saturation of the single diffractive cross section was 
described by Goulianos ||71|| using the term of a renormalized flux. 
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FIG. 3. The inelastic (production) cross section a^-Air calculated by DPMJET-II.5 as function of the laboratory collision 
energy (from 0.01 TeV up to 10^ TeV) compared to experimental data collected by Mielke et al. and to experimental data 
from the AGASA and Fly's eye |Q experiments as presented in corrected form in a recent paper by Block, Halzen and 
Stanev Given are also the calculated total and elastic cross sections. 
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FIG. 4. The inelastic (production) cross section ap^-N calculated by DPMJET-II.5 as function of the laboratory collision 
energy (from 0.01 TeV up to 1.E9 TeV). Given are also the calculated total and elastic cross sections. 



B. The Dual Parton Model for hadron production in hadron— nucleus and nucleus— nucleus collisions 



In the following we discuss new features in DPMJET-II.5, for a more detailed description of the model as applied 
in the code we refer to Refs. ||,||,|l|,||,|4|j3||3^,0,|9|j3O|||l. 



1. The Monte Carlo realization of the dual parton model DPMJET for hadron-nucleus and nucleus-nucleus collisions 



The model starts from the impulse approximation for the interacting nuclei - i.e. with a frozen discrete spatial 
distribution of nucleons sampled from standard density distributions | [7^ , |4^ . The primary interaction of the incident 
high-energy projectile proceeds via n elementary collisions between Hp ~ ua and rit = ub nucleons from the projectile 
(for incident hadrons rip — 1 ) and the target nuclei, resp, . Actual numbers n, Up and rit are sampled on the basis of 
Glauber's multiple scattering formalism using the Monte Carlo algorithm of Ref. [77| . 

The Glauber model, which is part of DPMJET-II allows to calculate the inelastic hadron-nucleus cross sections. 
The needed ingredients of this calculation are the nuclear geometry and the elementary hadron-nucleon scattering 
amplitude. 

The nuclear matter distribution is parametrized as follows: 

For most nuclei (untill DPMJET-II. 3/4 all nuclei) we parametrize the average nuclear radius (in fm) 



- - 1.12 V"-33 



(2.2) 



This parametrization is however not reliable for some light nuclei important for Cosmic Ray collisions. Therefore we 
replace (2.2) for these nuclei by the values given in Table 2 [|43| 



Table 2. Average nuclear radii used in DPMJET-II.5 |43|(in fm) 
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The nuclear density distribution is parametrized for most nuclei as the Fermi distribution 



Fir) = — — (2.3) 

However, for Na between 9 and 18 we use again more precise density distributions as given in Rcf. p3[ |. 
The elementary hadron-nucleus scattering amplitude is parametrized as follows: 

Energy dependent quantities enter the Glauber approach via the profile function of elastic hadron-nucleon scattering, 

7hN{b) = J (fqexp{iq- h)fhN{q), (2.4) 

i.e. the amplitude of elastic hadron-nucleon scattering in the impact parameter representation (with q denoting the 
lateral, i.e. two-dimensional momentum transfer). In their Monte Carlo realization of Glauber's approach Shmakov 
et al. 1?^ apply the high-energy approximation of the profile function, 

with parameters a*°^,a and p — 3fi//ijv(0) / 3//i7v(0) appropriate for the description of nucleus- nucleus interactions 
at energies of several GeV per nucleon. (This parametrization corresponds to a differential cross section da/dt ~ 
atot exp(a . t) with t ~ —q^-) 

However, the energy dependence of the elastic hadron-nucleon scattering amplitudes will influence the properties of 
hadron-nucleus and nucleus-nucleus scattering. In particular, the number of individual high-energy hadron-nucleon 
interactions {n,nA,nB) will increase with rising energy, hence the multiplicity will increase stronger than to be 
expected from the energy dependence of single hadron-hadron interactions. 

Guided by the data collected in Rcf. [|78[ , we apply the following parametrizations for the slope parameter a : 
a — 8.5 (1 -I- 0.0651ns) for nucleon-nucleon collisions above s = 50 GcV^, a = 6.2 (1 -I- 0.131ns) for nucleon-nucleon 
collisions below s = 50 GeV^ and a — 6.0 (1 -I- 0.0651ns) for tt- and K-nucleon collisions. We use for the ratio p 
of the real and imaginary part of the elastic scattering amplitude: p — —0.63 -|- 0.175 In for the energy region 
3.0 < ^/s < 50. and p — 0.1 in the energy region > 50 GeV in nucleon-nucleon scattering and p = 0.01 for tt- and 
K-nucleon scattering. 

The energy dependence of the total cross sections is described by the fits of the Particle Data Group ; at energies 
beyond the range of the actual parametrization of the pp cross section the one for pp is applied and at energies even 
higher we use the total cross-sections as calculated by the two-component Dual Parton Model for hadron-hadron 
collisions, see also Fig.|l|. 

The same information is also needed to construct the inelastic events and indeed for nuclear collisions each run of 
DPMJET starts with a calculation of the inelastic cross section. 

A new option (XSECNUC) has been added to DPMJET-II.5 to calculate a table of hadron-nucleus or nucleus- 
nucleus cross sections as function of the collision energy. The routine XSGLAU (this is a routine adopted from 
DTUNUC-II p^ , p5[ |) using the method according to |Q) calculates the following cross sections: 

(l)c'tot 

{2)<T^f^ A + B ^ A + B 

i^Kual-ei-s A + B^A + X, excluding (2) 

A + B^X + B, excluding (2) 
(5)<uaf.-ei-5 A + B^X, excluding (2), (3) and (4) 

('gW'4--B_ A-B_ A-B_ A-B_ rr^-^ _ a^~^ - (T^~^ 

^ / prod iiiel tot el quasi — el — Z quasi — el— A quasi — el — h 

In Fig. H and 1^ we present the cross sections (only crf^J^ , and cr^^^ ) as calculated by DPMJET-II.5 using 

the XSECNUC option for p-Air and Fe-N which are relevant for calculating the Cosmic Ray cascade as function of 
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the lab energy from 0.02 TeV to 10^ TeV. The quasi- elastic cross sections f^aS-ei-i only contained in these plots 
only in in af'^^ . The inelastic DPMJET p-Air cross sections cr'tzP are compared to the experimental data collected 
by Mielke et al. and to experimental data from the AGASAjT^ and Fly's eye experiments as presented in 
corrected form in a recent paper by Block, Halzen and Stanev [[75|. 



C. Implementation of new DPM diagrams for an improved description of baryon stopping in hadron— nucleus 

and nucleus— nucleus collisions 



1. Diquark fragmentation, the popcorn mechanism 



The fragmentation of diquarks is slightly more complicated than the fragmentation of a quark jet. As justified by 
Rossi and Veneziano ||80[| the baryon can be pictured as made out of three quarks bound together by three strings 
which join in a so called string junction point. In diagrams we can characterize the baryons 

(i) by the three quarks and the string junction or 

(ii) by a quark and a diquark, in this case the string junction always goes with the diquark. 

In all the diagrams discussed in this section we will either plot the quark and the diquark or, if the diquark breaks, 
the three quarks and the string junction. The quarks and the two quarks of a diquark are plotted as solid lines, the 
string junction is plotted as a dashed line. 

There are two possibilities for the first fragmentation step of a diquark, see FigJ|. Either we get in the first step 
a baryon, which contains both quarks of the diquark and the string junction or we get in the first step a meson 
containing only one of the two quarks and the baryon is produced in one of the following fragmentation step. This 
mechanism is well known, it is presented in the review on the Dual Parton Model and it was investigated for 
instance in |^,^ . This mechanism was implemented from the beginning in the BAMJET-fragmentation code 46 
used previously in DPMJET. This mechanism is also implemented under the name popcorn fragmentation in the Lund 
chain fragmentation model JETSET i3|,|3. 




FIG. 5. (a) Conventional diquark fragmentation, the baryon is produced in the first fragmentation step, (b) The popcorn 
mechanism: a meson is produced in the first fragmentation step, the baryon appears in the second or later fragmentation steps. 
We plot the lines for the two quarks of the diquark as solid lines, the string junction is plotted as a dashed line. 

What happens in the model with the popcorn mechanism compared to the model without can be most easily 
seen looking at the proton rapidity distribution in p-p collisions. The two maxima in the target and projectile 
fragmentation region of the proton rapidity distribution shift by about half a unity to the center, these maxima 
become wider and correspondingly the dip in the center is reduced. At the same time the Feynman x distributions of 
mesons get a component at larger Feynman x. The effects in hadron-nucleus and nucleus-nucleus collisions are quite 
similar. The popcorn mechanism is however not enough to explain the baryon stopping observed experimentally in 
hadron-nucleus and nucleus-nucleus collisions E8,E9], this will be discussed in more detail in Section III. 



2. New diquark breaking DPM diagrams in hadron-nucleus and nucleus-nucleus collisions 



New diquark breaking DPM-diagrams mainly of interest in hadron-nucleus and nucleus-nucleus collisions were 
proposed by Capella and Kopeliovich and investigated in detail by Capella and Collaborators [^;^ 85 1. Similar 
ideas were discussed by Vance and Gyulassy |86| . Capella and Kopeliovich [B0[ did discuss in detail their first diquark 
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breaking mechanism (see Fig.g, where this mechanism is characterized for nucleon-nucleon colhsions), in this case the 
valence-diquark breaks into the two quarks, the baryon is produced in the second or in later fragmentation steps. 




FIG. 6. The first C-K diquark breaking mechanism |5C|] plotted for a nucleon-nucleon collision. 

Looking only at the diagrams, one does not see any difference between the popcorn mechanism and the first Capella- 
Kopeliovich (C-K) diquark breaking mechanism. However, the mechanism of this first (C-K) diquark breaking 
mechanism differ in detail from the popcorn mechanism discussed above. Nevertheless, implementing the first C-K 
mechanism in DPMJET did not give any new feature of baryon stopping, which could not also be obtained from the 
popcorn mechanism. Therefore, we continue to use in DPMJET-1I.5 the popcorn mechanism instead of the C-K 
mechanism, there is no need to use both effects. 

More interesting is the second C-K mechanism, which was proposed but not discussed in detail in This 
mechanism was discussed in detail by Capella and Collaborators [pT| , p2|j8^ ] . In Figj^ we plot first the diquark- 
conserving diagram for a nucleon-nucleus collisions with two participants of the target nucleus. This is the traditional 
way for such a collision in the DPM. 




FIG. 7. The diquark-conserving diagram for a nucleon-nucleus collision with two participants of the target nucleus. 

In Fig.^ we plot the second C-K diquark-breaking diagram for the same collision. Now the second valence quark 
from the broken diquark in the first C-K diagram in Fig.^ is replaced by a Glauber sea quark from the nucleon 
projectile. Therefore, we will call the mechanism the Glauber sea quark mechanism of baryon stopping GSQBS. The 
probability of such a diquark splitting rises if the considered nucleon is involved in more than two interactions. The 
GSQBS has been implemented in DPMJET-II.5 and we will see in Section III, that in nucleon-nucleus collisions 
and nucleus-nucleus collisions we are able with this mechanism to fill the dip in the baryon rapidity distributions at 
central rapidity in agreement to the experimental data. As discussed already in detail in ||5^ , ^5|| this mechanism also 
contributes to increase the Hyperon production in nucleon-nucleus and nucleus-nucleus collisions. 
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FIG. 8. The Glauber sea quark mechanism of baryon stopping GSQBS for a nucleon-nucleus coUision with two participants 
of the target nucleus. This is the second C-K diquark-breaking mechanism [EH] . 



3. The Casado diagram 

A new diagram was also introduced by Casado ||53|1 . The diagram is plotted again for a nucleon-nucleus collision 
with two participants of the target nucleus in Fig. M. Here the fragmenting diquark contains a Glauber sea quark and 
hadronizes like a valence diquark producing baryons mainly in the fragmentation region. However, the flavor content 
of the baryon is changed. This diagram gives another contribution to hyperon production in nuclear collisions. We 
have implemented the Casado diagram in DPMJET~II.5 and it is used in addition to the GSQBS diagram for all 
nuclear collisions compared in Section III to data. 




FIG. 9. The diagram introduced by Casado for a nucleon-nucleus collision with two participants of the target nucleus. 
The diquark contains one Glauber sea quark. 

4- Two extrapolations of DPMJET-II.5 to Cosmic Ray energies. 




FIG. 10. Standard DPM diagram for a nucleon-nucleon interaction with one additional soft secondary interaction induced 
by the unitarzation procedure. There is one valence-valence and one sea-sea interaction, each represented by a pair of chains. 
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FIG. 11. New DPM diagram for a nucleon-nucleon interaction with one additional soft secondary interaction induced by the 
unitarzation procedure. The diquark is split and an unitarity sea quark is used similar to Fig. ^ to shift the baryon in one of 
the chains. We call this the unitary sea quark mechanism for baryon stopping USQBS. 

At high energies we have multiple collisions even in hadron-hadron collisions due to the unitarization procedure. 
We call the sea quarks at the ends of the additional chains in this case unitary sea quarks. The Glauber sea quarks 
are needed in nuclear collisions already at rather low energies, for instance at the energies of heavy ion collisions at 
the CERN-SPS. In contrast to this, unitary sea quarks appear in significant numbers in hadron-hadron and nuclear 
collisions only at rather high energies, for instance at the energies of the CERN-SPS collider or the TEVATRON 
collider, they are important in the Cosmic Ray energy region. 

With the unitary sea quarks at the ends of the chains from the secondary collisions we obtain a new mechanism for 
baryon stopping, which will become effective at very high energies. 



In Fig, 10 we plot the standard DPM diagram for a nucleon-nucleon interaction with two soft interactions induced 
by the unitarization procedure. There is one valence-valence and one sea-sea interaction, each represented by a pair 
of chains. In analogy to Fig. ^ we construct from this the new diagram for baryon stopping in Fig. The diquark 
is split and an unitary sea quark is used to have the baryon only in the second or later fragmentation steps in one of 
the chains. We call this the unitary sea quark mechanism for baryon stopping USQBS. Also here the probability for 
such a diquark splitting rises if there are more than 2 interactions of the hadrons involved. Obviously, this mechanism 
leads to Feynman x distributions of baryons in p-p collisions becoming softer and Feynman x distributions of mesons 
becoming harder than without the USQBS mechanism. 

In contrast to the GSQBS mechanism, which leads already to effects in nuclear collisions at the energy of the CERN- 
SPS, we have at present no data available to prove that this USQBS mechanism is a valid extension of the DPM. There 
is no data for baryon stopping in proton-proton or antiproton-proton collisions at collider energies. Unfortunately, 
the fragmentation region at large Feynman x has not been investigated experimentally with enough detail at any of 
the hadron-hadron colliders. However, if the GSQBS mechanism is the correct mechanism responsable for the baryon 
stopping effects found in nuclear collisions, then also the USQBS mechanism should modify the collisions at collider 
energies and beyond. In Section IV, where we present the properties of DPMJET-II.5 at Cosmic Ray energies we 
present always two extrapolations of the model. In one version (shortly characterized as version 55 in the plots) we 
have both the GSQBS and USQBS mechanism, in a second version (characterized in the plots as version 50) we have 
only the GSQBS mechanism. Only collider experiments on baryon stopping or Cosmic Ray observations can prove, 
which version is the better one. From theoretical prejudices however there would be the claim, that the version with 
both the GSQBS and USQBS mechanisms is to be prefered. 



5. New parameters connected with the diquark breaking diagrams 

For each of the new diquark breaking diagrams described in this Section we have to introduce a new parameter. 
These parameters give the probability for the diquark breaking mechanisms to occur, given a suitable sea quark is 
available and given that the diquark breaking mechanism is kinematically allowed. For an original diquark-quark 
chain of small invariant mass, which originally just fragments into two hadrons, the diquark breaking is often not 
allowed at small energies. 

The optimum values of the new parameters are determined by comparing DPMJET-II.5 with experimental data 
on net-baryon distibutions. We get roughly for all of these probabilities values around 0.5. The definitions of these 
parameters and their default values in DPMJET-II.5 are given in the companion paper Q. 
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D. Realization of the intranuclear cascade process 



Nothing has been changed in DPMJET-II.5 against DPMJET-II.3/4 with respect to the formation zone intranuclear 
cascade, nuclear evaporation and fragmentation, 7 deexcitation of the residual nuclei and the Cronin effect. All of 
this was described in detail in ||^,|2^,^ . 



E. Production of strange particles 



Studies of strangeness production within this model were given in [ p5| , p7y3§| ]. Enhanced generation of strange 
particles, in particular of strange antibaryons, has been proposed as a signal for the formation of quark gluon plasma 
in dense hadronic matter Recent data from experiments at the CERN SPS have already been interpreted 

within this scheme. However, we find it worthwhile to pursue the study of tested conventional models without QGP 
formation like the DPM before drawing final conclusions. The DPM is an independent string model. Since the 
individual strings are universal building blocks of the model, the ratio of produced strange particles over non-strange 
ones will be approximately the same in all reactions. However, since some strings contain sea quarks at one or 
both ends and since strange quarks are present in the proton sea, it is clear, that, by increasing the number of 
those strings, the ratio of strange over non-strange particles will increase. This will be the case for instance, when 
increasing the centrality in a nucleus-nucleus collision. It is obvious, that the numerical importance of the effect 
will depend on the assumed fraction of strange over non-strange quarks in the proton sea. The rather extreme case 
leading to a maximum increase of strangeness would be to assume a SU(3) symmetric sea (equal numbers of u, d 
and s fiavors). We express the amount of SU(3) symmetry of the sea chain ends by our parameter s"*^" defined as 
gsea — 2 < Ss >/(< Us > + < ds >) whcrc the < qs > give the average numbers of sea quarks at the sea chain ends. 
Usually, DPMJET uses s"""- = 0.5. 

There are a number of effects, which change the number of strange hadrons, especially strange and multistrange 
baryons and antibaryons in the model: 

(i) The presence of sea qq and qq diquarks at sea-chain ends | ]89| , p0| . 

(ii) Thc popcorn effect. 

(iii) The diquark breaking diagrams discussed in Section H.C, see also 52 8£ 



(iv) The secondary interactions of co-moving produced particles ||9l|,|92|, sec also Section H.G. 

(v) Effects of string fusion ||93| , |9^ and percolation of strings l95| , |9q| . 

It is certainly beyond the scope of the present paper, to discuss all of these effects in detail. The effects (i) to (iv) 
are contained in DPMJET-II.5. 



F. Single diffractive cross sections in hadron— nucleus and nucleus— nucleus collisions 

Let us first summarize the sampling of single diffractive events in in DPMJET in hadron-hadron collisions. 

Single diffraction within the Dual Parton Model was studied in detail and compared to experimental data in p6| , [76| . 
Single diffraction dissociation is represented by a triple-Pomcron exchange (high mass single diffraction) and a low 
mass component (low mass single diffraction) JTst . 

Diffractive processes characterized by the excitation of an initial hadron to intermediate resonances (low mass 
diffractive interactions) are introduced via a two channel eikonal formalism. 

High mass single diffractive events are sampled cutting the triple-Pomeron graph . The excited system consists of 
two chains stretched between the constituents of the excited hadron and a qq-pair of the Pomeron emitted from the 
hadron at the upper vertex. 

Corresponding to M|, = xds and to the known 1/M^-behaviour of the diffractive differential cross section 
d^asd/dtdM"^ we sample the sum of the momentum fractions of the qq-pair xd — Xq + Xq from a l/xc-distribution. 
Assuming that inelastic diffraction dissociation puts virtual hadronic states on the mass shell we limit this selection 
by the coherence condition. 

Ml/s<{mpRy^ with R ~ m-'^ (2.6) 
where ^/s is the cm. energy of the collision and R is the interaction radius. 

An excited low mass single diffractive event is represented by a chain connected to the valence partons of one hadron. 
The invariant mass Md is fixed in the same way as above. 

We turn to single diffractive cross sections in hadron-nucleus and nucleus-nucleus collisions. 
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In DPMJET the simple approximation is used: single diffraction in nuclear collisions occurs only in Glauber 
collisions with only one participant of projectile as well as target {Na = Nb = N = 1). In such collisions also 
the total number of Glauber collisions is equal to one. In the first step of the Na — Nb = N ~ 1 collision we 
have just a single hadron-hadron collision, single diffraction occurs, if this single collision is a diffractive collision. 
In the second step of the single diffractive hadron-nucleus or nucleus-nucleus collision we have the formation zone 
cascade in the participating nuclei followed by nuclear evaporation. There might be single diffractive collisions also 
in Glauber collisions with N larger than one. In such cases it is rather unlikely, that all collisions will be diffractive. 
If we have one single diffractive collision together with one or several nondiffractive collisions, the event just looks 
like a nondiffractive one. We neglect in DPMJET the small fraction of diffractive events, where all of these Glauber 
collisions are diffractive. 

DPMJET allows in hadron-nucleus and nucleus-nucleus collisions like in hadron-hadron collisions to chose all 
collisions or only nondiffractive or only diffractive collisions. DPMJET always calculates the fraction of diffractive 
collisions to the inelastic cross section. 

In Fig.|l^ we plot the single diffractive cross sections in proton-nucleus collisions as function of the nuclear mass 
number A. The single diffractive cross sections given are always for diffraction on the target as well as the projectile 
side. As to be expected we find cr^^)^ ~ A^^^. Single diffractive nucleus-nucleus cross sections (Jg^^ are found to be 
rather similar to the h-Max(A,B) diffractive cross sections. In Table 3 we give only some examples of single diffractive 
nucleus-nucleus cross sections obtained from DPMJET-II.5. 



sigpadiff200 



100 - 



SD 
mb 




10 100 

Target mass number 

FIG. 12. Single diffractive cross sections calculated with DPMJET-II.5 in proton-nucleus collisions. The single diffractive 
cross sections given are for diffraction on the target as well as the projectile side. 

Table 3. Diffractive cross sections in nucleus-nucleus collisions (in nib). The single diffractive cross sections given 
are for diffraction on the target as well as the projectile side. 



Reaction 




S~S 


50.5 


S-Ag 


60.9 


S-Au 


67.6 


Pb-Pb 


72.8 
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G. Final state interactions of co— moving secondaries in nuclear collisions 



Analyzing the rapidity distributions of produced A and A perticles in central heavy ion collisions at CERN-SPS 
energies within the Dual Parton Model Capella et al noted the need for secondary interactions of co-moving 

secondaries to understand the data. They only introduced the following secondary interactions 

TT + N ^K + K (2.7) 

and 



TT + N ^ K + K (2.8) 



A reasonable cross section for these reactions of a ~ 1.5 mb was needed in 92| to understand the data. All of 
this was done in |9^:^^ using analytical methods, not a Monte Carlo event generator. Final state interactions of 
co-moving secondaries was also used in [0,|8| to explain the J/tp supression in Pb-Pb collisions. 

The method of Capella et al ||l|,||| was implemented in DPMJET aheady in 1995. In DPMJET the method is used 
to modify the Monte Carlo events and we use only a cross section of cr « 1.0 mb. First results reproducing essentially 
the results of [ ^l]j9^ are described shortly in and in more detail in a unpublished code write-up Q. We would 
like to stress, this method as implemented with only the two reactions given above, can only be considered as a first 
preliminary step. A better method of final state interactions should be implemented as a Monte Carlo method from 
the beginning and more types of secondary interactions should be taken into account including the reverse reactions. 
With the present method we can not expect to obtain reasonable results in situations with much higher secondary 
particle densities than in the heavy ion collisions at the CERN-SPS energy. Therefore, it is not recommended to use 
the secondary interaction option of DPMJET at RHIC or CERN-LHC energies. 



H. Central nucleus— nucleus collisions in DPMJET 



The sampling of a Glauber event according to the method of proceeds as follows: First the impact parameter 
B of the hadron-nucleus or nucleus-nucleus collision is sampled from a distribution determined at the beginning by 
a Monte Carlo method. In a second step the number of interacting projectile and target nucleons Na — Np and 
Nb = Nt and the total number of Glauber collisions N are sampled. 

We might define in DPMJET a central collision by different methods: 

(i) Allowing only impact parameters B < Bthr- 

(ii) Demanding, that a certain minimum fraction of the projectile or target nucleons takes part in the collision for 
instance Na = Np > aA, where a is a number not far from unity. 

Experiments investigating central nucleus-nucleus collisions are usually not able to determine reliably the impact 
parameters of the collisions. They might be able to determine the number of interacting projectile nucleons, for 
instance using a zero-degree calorimeter. Most often one finds however the central collision characterized by giving 
which fraction Fcentrai of all events are contained in the sample of most central collisions considered. 

In order to keep contact to the Fcentrai definition we proceed in DPMJET in an iterative way. DPMJET is run 
rejecting all events with Na = Np < aA for a given a. At the same time we keep track which fraction of the events 
fa is accepted in the run. a is iterated untill fa « Fcentrai- 



I. Restrictions of DPMJET-II 



DPMJET-II.5 works for hadron-hadron and hadron-nucleus collisions above 50 MeV kinetic energy, it has been 
compared to data on average multiplicities down to about 1 GeV and to data on Feynman x distibutions down to 
3-4 GeV see The HADRIN model used by DPMJET for hadron-hadron collisions below 5 GeV was compared to 
data down to several hundred MeV ||9£,10C]. DPMJET-II.5 works for nucleus-nucleus interactions above 5 GeV per 
nucleon, but was never compared to experimental data below 60 GeV. 

DPMJET-II.5 is able to run up to energies of approximately 10^^ eV in the lab system (or -/s — 2000 TeV). 
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III. COMPARING DPMJET-II.5 TO DATA 



A. Comparing to data in hadron— hadron collisions 



Before a model like DPMJET-II.5 can be used for anything, we have to demonstrate, that the model describes well 
enough the available data up to the highest energies. The parameters of the model as given in Table 1 are different 
from the parameters used in the previous version of the model, therefore, this agreement to the data is not trivial, 
even if the previous DPMJET version was quite successful in this respect. 

Table 4. Comparison of avera ge m ultiplicities of produced hadrons in proton- proton collisions at 200 GeV. The 
experimental data are from Ref. |101|. 



Particle 


DPMJET-II.5 


Exp. |101] 


rich 


7.69 


7.69 ± 0.06 


ri- 


2.85 


2.85 ± 0.03 


P 


1.24 


1.34 ± 0.15 


n 


0.68 


0.61 ± 0.30 


7:+ 


3.25 


3.22 ± 0.12 




2.58 


2.62 ± 0.06 




3.55 


3.34 ± 0.24 


K+ 


0.30 


0.28 ± 0.06 


R- 


0.19 


0.18 ± 0.05 


Kl 


0.23 


0.17 ± 0.01 


A 


0.14 


0.096 ± 0.01 


A 


0.0178 


0.0136 ± 0.004 


P° 


0.48 


0.33 ± 0.06 


P 


0.048 


0.05 ± 0.02 



15 



We start with looking at data, which prove, that the average multipUcity and the particle content are well described. 
In T able 4 we compare to the very well known multiplicities of secondary hadrons in 200 GeV p-p collisions from Ref. 
1 101 1 . In Table 5 w e com pare to multiplicities of secondary hadrons and hadron resonances in -^s = 27.5 GcV p-p 
collisions from Ref. |l02| . 

Table 5. Comparison of average multiplicities of produced hadrons and hadron resonances in proton- proton 
collisions at — 27.5 GeV. The experimental data are from Ref. |102|. 



Particle 


DPMJET-II.5 


Exp. pO^ 


P 


1.28 


1.20 ± 0.097 ± 0.022 


7T+ 


3.52 


4.10 ± 0.11 ± 0.15 




2.87 


3.34 ± 0.08 ± 0.12 


ttO 


3.88 


3.87 ± 0.12 ± 0.16 


K+ 


0.34 


0.331 ± 0.016 ± 0.007 


R- 


0.24 


0.224 ± 0.011 ± 0.004 


K% 


0.27 


0.232 ± 0.011 ± 0.004 


A 


0.156 


0.125 ± 0.008 ± 0.008 


A 


0.026 


0.020 ± 0.004 ± 0.0008 


S+ 


0.050 


0.048 ± 0.015 ± 0.004 




0.016 


0.0128 ± 0.0061 ± 0.0032 


A++ 


0.24 


0.218 ± 0.0031 ± 0.013 


AO 


0.23 


0.141 ± 0.0098 ± 0.0089 


A++ 


0.005 


0.013 ± 0.0049 ± 0.0049 


A" 


0.006 


0.0336 ± 0.008 ± 0.0006 


P° 


0.54 


0.385 ± 0.018 ± 0.038 


P+ 


0.54 


0.552 ± 0.083 ± 0.046 




0.43 


0.355 ± 0.058 ± 0.033 




0.50 


0.390 ± 0.024 ± 0.002 


V 


0.55 


0.30 ± 0.02 ± 0.054 


K*+ 


0.18 


0.132 ± 0.016 ± 0.002 


K*- 


0.11 


0.088 ± 0.012 ± 0.001 




0.15 


0.119 ± 0.021 ± 0.002 




0.11 


0.0903 ± 0.016 ± 0.001 




0.029 


0.019 ± 0.0018 ± 0.054 


P 


0.066 


0.063 ± 0.002 ± 0.001 
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FIG. 13. Comparison of Feynman-x distributions of positively charged hadrons produced in proton-proton collisions at 250 
GeV. The experimental data are from the EHS-NA22 Collaboration [L03|. 
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FIG. 14. Comparison of Feynman-x distributions of vr 



experimental data are from the EHS-NA22 Collaboration [ Wc ] 
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mesons produced in proton-proton collisions at 250 GeV. The 
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FIG. 15. xiab cross sections in p-p collisions. We compare data at different energies 104,102,10E] shown as symbols with 



DPMJET-II.5 at piab = 200 GeV/c shown as sohd triangles. The data are also compared to the SIBYLL |10^ and QGSJET 
[L07| models shown as histograms. 
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FIG. 16 . Ener gy fraction xiab carried by the leading proton. The data are photoproduction and DIS measurements at - 



200 GeV [106,109] shown as symbols compar ed to DPMJET-II.5 for p-p collisions at ^/s — 200 GeV shown as open triangles. 
The data are also compared to the SIBYLL [106 and QGSJET [107[ models shown as histograms. 
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In Fig. and |Ij we compare the model with NA-22 data on the Feynman-x distribution of positively charged 
hadrons and of produced in 250 GeV pp collisions. In Fig. ^ we observe in particular (compared to previous 
DPMJET versions) a better agreement of the model predictions with the forward production of protons. 

and|l|. The Figs. 
The leading 
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Distributions of leading protons are com pared more directly to data in the following two Figs, 
are adopted from talks of Engel |11C,111|. But of course, we present now the DPMJET-II.5 results, 
particle production is very important for the Cosmic Ray cascade simulation. 

In Fig. [l^ we compare xiab cross sections in p-p collisions measured at different energies [104,102 105 1 with 
DPMJET-II.5 results at piab = 200 GeV/c. The DPMJET cross sections change only very slowly with energy. The 



obtained agreement with DPMJET-II.5 to the d ata is much be tter than with former DPMJET versions see [110|. 
The data are also compared to the SIBYLL [|0§ and QGSJET |l07) models. 

In Fig. ^ we compare the energy fraction xiab carried by the l eading p roton. The data are photoproduction and 
DIS measurements from the HERA Colhder at Vi = 200 GeV p8|jl09| . We compare to DPMJET~II.5 for p-p 
collisions at y/s — 200 GeV. The forward production of leading protons is not expected to depend strongly on the 
reac tion channel. It is found again, that DPMJET-II.5 agrees much better to the data than older DPMJET versions, 
see 111 . This comparison demonstrates, that the leading proton distribution at y/s — 200 GeV is still rather flat like 
at the much lower energies in Fig. |l|. At = 200 GeV both versions of DPMJET-II.5 with GSQBS and USQBS 
and with only GSQBS, (discussed in Section II) lead still to very similar xiab distributions of secondary protons. 
Unfortunately therefore, these data cannot be used to discriminate between the two DPMJET-II.5 versions. The 



data are also compared to the SIBYLL flO^ and QGSJET ||107[ models 
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FIG. 17. Com parison with transverse momentum cross sections at ^/s = 0.2 and -.9 TeV with collider data from the UA- 
Experiment [L12|. 
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FIG. 18. Co mpa rison witli transverse momentum cross sections at ^/s = 0.63 TeV with collider data from the UA-1 MIMI 
Collaboration [113]. The experimental data are represented by the parametrization given by the Experiment. 
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FlC. 19. Average transverse momenta of charged secondaries according to the two versions of DPMJET-II.5 produced in pp 
colhsions in the full pseudorapidity interval as function of the center of mass energy ^/s. We compare to data collected by the 
UA-1 Collaboration [112] which refer to the central pseudorapidity range and are therefore expected to be slightly larger than 
in the full range. 



Hadron transverse momentum distributions were measured by the UA-1 Collaboration |112|. In Fig.y/7|we compare 
the distributions at y/s = 200 GeV and 900 GeV with DPMJET-II.5. The transverse momentum distribution up to 
larger p± was determined by the UA-l-MIMI Collaboration ||lT|l. In Fig.|l| we compare the model results with the 
parametrization of the data given by this experiment. 
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In Fig.[T^ we compare average transverse momenta as obtained from the two versions discussed in Section II of 
DPMJET-II.5 as function of the cms energy ^/s with data collected by the UA-1 Collaboration. This plot gives at 
the same time the DPMJET-IL5 predictions for the average transverse momenta up to y/s = 1000 TeV. We find, that 
the average transverse momenta in both versions of the model agree rather well to each other. Average transverse 
momenta depend on the longitudinal variables, they differ if one chooses different pseudorapidity ranges, the data 
given for the central pseudorapidity range are expected to be slightly larger than the calculation given for the full 
range. 

With rising energy the fraction of strange hadrons rises slowly. In Fig. ^ th e K/ TT ratio according to DPMJET-II.5 
for pp or pp collisions is compared to data from the E735 Collaboration |114|. 
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FIG. 20. K/tt ratios in pp or pp collisions as function o f the cms energy ^/s. The DPMJET-II.5 calculation is compared with 
data collected from the E735 Collaboration at Fermilab [114|. 
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B. Comparing to data in hadron— nucleus and nucleus— nucleus collisions 

We now turn to collisions with nuclei. 

In Table 6 we compare average multiplicities of negatively charged hadrons calculated with DPMJET-II.5 in 
hadron-hadron, minimum bias hadron-nucleus and central nucleus-nucleus collisions with experimental data. 

Table 6. Comparison of average multiplicities of produced negatively charged hadrons in proton-proton, minimum 
bias proton-nucleus and central nucleus-nucleus collisions at 200 GeV. 



Collision 


DPMJET-II.5 


Exp. 


Reference 


P-P 


2.85 


2.85 ± 0.03 


1 


101| 


p-S 


5.10 


5.0 ±0.2 




m 


p-Ar 


5.30 


5.39 ±0.17 






p-Ag 


6.18 


6.2 ±0.2 


1 


117| 


p-Xe 


6.43 


6.84 ±0.13 






p-Au 


6.81 


7.0 ±0.4 


1 


117| 


p-Au 


6.81 


7.3 ±0.3 






S-S central 


103 


98 ±3 


|9 




S-Ag central 


174 


186 ±11 


|9 




S-Au central 


202 


225 ±12 


H 





Table 7. Comparison of average multiplicities of produced strange hadrons in central S-S and S~Ag collisions at 
200 GeV. For DPMJET-II.5 we give the results for the models without and with secondary interactions of com overs 
and it is assumed that I]° and S*^ have decayed. The experimental data are from the NA35 Collaboration ]lig| ]. 



Reaction 


Particle 


DPMJET-II.5 


DPMJET-II.5 with sec.int. 


Exp. [119| 


S-S 


A 


5.7 


7.3 


9.4 ± 1.0 


S-S 


A 


1.13 


1.22 


2.6 ± 0.3 


S-S 




8.7 


9.1 


10.5 ± 1.7 


S-S 


K+ 


10.6 


11.4 


12.5 ± 0.4 


S-S 


K- 


6.8 


6.9 


6.9 ± 0.4 


S-Ag 


A 


9.6 


13.1 


15.2 ±1.2 


S-Ag 


A 


1.8 


2.0 


2.6 ± 0.3 


S-Ag 




14.6 


15.5 


15.5 ±1.5 


S-Ag 


K+ 


17.8 
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In Tables 7 and 8 we give the average multiplicities of produced strange hadrons in central S-S, S-Ag and Pb-Pb 
collisions at 200 GeV (158 GeV for Pb-Pb). We give the DPMJET-II.5 results for the models without and with 
secondary interactions and it is assumed, that YP and E° have decayed. The experimental data for central S-S and 
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S-Ag collisions are from the NA35 Collaboration |119|. There is some doubt, that the simple model for secondary 
interactions in DPMJET-II.5 gives correct results in a situation with such high produced hadron density like in the 
central Pb-Pb collisions. 

In Table 8 we give the results for central Pb-Pb collisions also for the RHIC and LHC energies and we include 
charged multiplicities and and central plateau heights in the Table. At the LHC DPMJET-II.5 gives a slightly lower 
plateau then DPMJET-II.2 (for which the plateau was given in |12C|). The central plateaus at RHIC and LHC 
energies given in Table 8 agree however very well with the ones given by Capella in . 

Table 8. Comparison of average multiplicities of produced strange hadrons in central Pb-Pb collisions at 158 GeV 
and at the energies of RHIC and the LHC. For DPMJET-II.5 we give only at 158 GeV the results for the models 
without and with secondary interactions of comovers and it is assumed that E*^ and have decayed. 
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In Fig. ^.a the comparison is with the rapidity distribution of negatively charged hadrons in p-p, p-Ar and p-Xe 
collisions at 200 GeV and in Fig. ^b we compare to rapidity distributions of negatively charged hadrons in central 
S-S and S-Ag collisions. 
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FIG. 21. (a) Negatively charged particle rapidity distribution for p-p, p-Ar and p-Xe interactions. The DPMJET-II.S 
results are compared with data Ml . (b) Rapidity distribution of negatively charged hadrons in central S-S ans S-Ag collisions. 
The results of DPMJET-II.5 are compared with data from the NA-35 Collaboration 
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In Fig. ^ we compare rapidity distributions of 7r+, tt", and K~ mesons produced in central (4%) Pb- 
Pb collisions a t 15 8 AGeV/c. The experimental data, mainly in the central rapidity region are from the NA-44 
Collaboration [[l2l[. The agreement of model and data is encouraging, however not perfect. 
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FIG. 22. (a) Rapidity distributions of tt"*" and mesons produced in cen tral Pb-Pb collisions. The DPMJET-ILS 
distributions are compared with experimental data from the NA-44 Collaboration |121|. (b) Rapidity distributions of n~ and 
K~ mesons produced in c entra l Pb-Pb collisions. The DPMJET-II.5 distributions are compared with experimental data from 
the NA-44 Collaboration |l2l| . 
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C. Baryon stopping 



We present first the DPMJET predictions at piab — 200 GeV/c for net baryon rapidity distributions in the original 
model without the new GSQBS and USQBS diagrams. In Fig. |2^.a we present the leading net proton (p — p) rapidity 
distribution dNp/dy — dNp/dy in p-p, p-S, central S-S and central Pb-Pb collisions. In Fig. 23, b we present the 
net A (A — A) rapidity distribution dN\/dy — dNj^/dy in p-p, p-S, central S-S and central Pb-Pb collisions. In 
p-p collisions, which at the given collision energy are hardly modified by the new GSQBS and USQBS diagrams, we 
observe a dip at central rapidity. This dip is also present in DPMJET without the new GSQBS and USQBS diagrams 
in p-S and central A-A collisions. This disagrees to the data presented in the next Figures. 
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FIG. 23. (a) Net proton (p — p) rapidity distribution in p-p, p-S, central S-S and central Pb-Pb. collisions. Calculated 
with DPMJET-II.5 without the new diagrams modifying baryon stopping, (b) Net A (A — A) rapidity distribution in p-p, 
p-S, central S-S and central Pb-Pb. collisions. Calculated with DPMJET-II.5 without the new GSQBS and USQBS diagrams 
modifying baryon stopping. 
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In Fig. |2J.a and we compare the net-proton distributions according to the full DPMJET-II.5 model with data 
in p-S and p-Au collisions p9|] . Now the dips at central rapidity are filled in the model, we observe like in the data 
at central rapidity a flat net-proton rapidity distribution. 
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FIG. 24. (a) Net proton (p — p) rapidity distribution in p-S collisions. The DPMJET-II.5 results are compared with data 
49[. (b) Net proton (p — p) rapidity distribution in p-Au collisions. The DPMJET-II.5 results are compared with data P9|. 
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In Fig. ^ we compare the full DPMJET-II.5 model in the two versions without and with secondary interactions 
with data on net -proton production in central S-S collisions. Also here the dip at central rapidity in the model 
has disappeared (compare to Fig.^^, however, the agreement to the data is not perfect. There is a significant 
disagreement in the fragmentation regions of the two nuclei. The reason for this is in the spectator evaporation 
protons included in DPMJET, but apparently not in the data. 
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FIG. 25. Net proton (p — p) rapidity distribution in central S-S collisions. The DPMJET-II.5 results are compared with 
data [^. (a) DPMJET-II.5 without secondary interactions, (b) DPMJET-II.5 with secondary interactions. 
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In Fig. we compare the full model to data on net-A production in p-S collisions. In Fig, 26 the dip at central 
rapidity has disappeared in the model. In Fig. ^a and .b we compare the full model to data on net-A production 
central S-S collisions for the two models with and without secondary interactions. In Fig, 27, a and .b again the central 
dips in the models are largely reduced but the agreement to the NA35-data is rather unsatisfactory in the model 
without secondary interactions. Also the agreement of the model without secondary interactions to the total number 
of A's produced according to NA35 is rather unsat isfac tory, see Table 7. However, we find no disagreement to the 
total-A rapidity distribution as measured by Na36 [122|. 

We might conclude from the comparison of net-p and net-A production, that the secondary interaction is indeed 
needed to understand the data, further improvements to the secondary interactions might be useful. 
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IV. PROPERTIES OF THE MODEL IN THE HIGHEST ENERGY REGION 



In all plots in this Section we present the DPMJET-II.5 results for the two versions of the model discussed already 
in Section II. C: 

(i) With only the new diagram for baryon stopping GSQBS, this version is characterized in the plots by 50, this 
version corresponds in its high energy behaviour approximately to the previous version of the model DPMJET-II.4. 

(ii) With both new diagrams for baryon stopping GSQBS and USQBS, this version is characterized in the plots by 
55, the hight energy behaviour of this version is new, at high energies the baryons carry 5 to 10 percent less energy 
than in version 50, correspondingly the mesons carry more energy than in version 50. 

Fig.^ presents the rise of the total charged multiplicity in p-p and p-N collisions with the cms energy ^/s. 
The corresponding rise of the average transverse momenta was already presented in Figjis] 
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Following for instance the basic discussion of |123|, we introduce a variable xiab similarly to Feynman-xp . but this 
time in the lab-frame : 



Xlab 



Eo 



(4.1) 



Ei is the lab-energy of a secondary particle i and Eq is the lab-energy of the projectile in a h-nucleus collision. 
We introduce xiab distributions F(xiab) ■ 
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FIG. 29. Spectrum weighted moments for pion production in p-p and p-N collisions as function of the (nucleon-nucleon) 
cms energy -^s. 



Fiixiab) = Xiab 
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(4.2) 



We note that the Feynman-Xi;- distribution at positive xp in the projectile fragmentation region is a very good 
approximation to the xiab distribution. 
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The cosmic ray spectrum-weighted moments in p-A colhsions are now defined as moments of the F{xiab) ■ 



n 



(4.3) 



Here —7 ~ -1.7 is the power of the integral cosmic ray energy spectrum and A represents the target nucleus. 
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FIG. 30. Spectrum weighted moments for Kaon production in p-p and p-N collisions as function of the (nucleon-nucleon) 
cms energy -^s. 



The spectrum-weighted moments for nucleon-air collisions, as discussed in [ 123 |, determine the uncorrelated fluxes 
of energetic particles in the atmosphere. 

We also introduce the energy fraction Kf~'^ : 



Kf ^ ^ I Ff ^{xiab)dxiab 
Jo 



(4.4) 



As for xiab, the upper limit for X is 1 in h-nucleus collisions. 
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In Figs.|2^ we present the spectrum weighted moments summed over pions of both charges in pp and p-N coUisions 
as function of the cms energy ^/s per nucleon. In Figs.^ the moments are given for charged Kaon production also in 
pp and p~N coUisions. 
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In Fig.^ we present again for pp anf p-N collisions the energy fractions K for net baryons B — B (baryon minus 
antibaryon), B (baryon) and charged pion production. The energy fraction Kg_g is always smaller than the energy 
fraction Kg. The difference between both is the energy fraction going into antibaryons Kj^ (not given separately) 
which is equal to the energy fraction carried by the baryons which are produced in baryon-antibaryon pairs. 

We find in DPMJET-II.5 all average values characterizing hadron production: the cross sections (Fig.[l|), the average 
transverse momenta (Fig.|l^) the charged multiplicities (Fig.^8|), and the moments in Figs. |3y and ^ to change 
smoothly with energy in most cases just like the logarithm of the energy. 

So far, we have not found any experimental data to favour either the versions of DPMJET-II.5 (i or 50) with only 
GSQBS or (ii or 55) with GSQBS and USQBS. Clearly, the version prefered by theoretical prejudices is (ii or 55). In 
this version we have a better Feynman scaling of meson distibutions and spectrum weighted moments and we have a 
faster decrease with the collision energy of the energy fractions into secondary baryons. 
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V. SUMMARY 



In the present paper we present new features of the DPMJET-II.5 model and the properties of the model in the 
Cosmic Ray energy region. The most important of the new features are the new GSQBS and USQBS diagrams, 
which modify significantly the baryon stopping in the model. Using as well the GSQBS and USQBS options we find 
significant differences in the high energy behaviour of the new model against former versions of DPMJET. 

We report also about comparisons of DPMJET-II.5 with accelerator and collider experiments. We find especially 
an improved agreement of the Feynman x distributions of leading protons to data and a reasonable agreement to 
experimental data on leading baryon stopping in nuclear collisions. Many of the published comparisons of former 
versions of the model to other aspects of experimental data come out quite similar with DPMJET-II.5, we mention 
especially the properties of the model at low energies [^ ^9||30[| and charm production |^ . 

There are certainly more tests of the model needed and new features to be incorporated for hadron-nucleus and 
nucleus-nucleus collisions, for strangeness production and for secondary interactions of the produced particles. DP- 
MJET will continue to evolve, there are already plans for a major new update to DPMJET-III. 
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